Exploring binding energy and separation energy dependences of HBT strength 
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Hanbury Brown-Twiss (HBT) results of the nucleon-nucleon correlation function have been pre- 
sented for the nuclear reactions with neutron-rich projectiles (Be isotopes) using an event-generator, 
the Isospin-Dependent Quantum Molecular Dynamics model. We explore that the relationship be- 
tween the binding energy per nucleon of the projectiles and the strength of the neutron-proton HBT 
at small relative momentum. Moreover, we reveal the relationship between the single neutron sepa- 
ration energy and the strength of the halo neutron-proton HBT. Results show that neutron-proton 
HBT results are sensitive to binding energy or separation energy. 
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I. INTRODUCTION 



Hanbury Brown-Twiss (HBT) technique was presented 
for the astrophysical measurements a few decades ago, 
and it reveals information about the angular diameters 
of distant stars 0. More recently, it has been widely 
used in others fields, for instance, the analogous correla- 
tions in semiconductors and in free space aiming at the 
fermionic statistics of electrons It has also be- 

come an important tool in high energy region since it can 
be utilized to measure the evolving geometry of the inter- 
action zone while being applied to the studies to search 
for a possible quark-gluon plasma and study the proper- 
ties of the predicted new state of matter 5] . In the past 
several years it has also significant theoretical develop- 
ment and widespread application in subatomic physics 
H, 0, H, 0, mi US ■ The emission time and source size 
in the nuclear reaction can be extracted by the nuclear 
HBT technique, which is one kind of intensity interfer- 
ometry. In the applications of experimental and theo- 
retical heavy ion reactions at intermediate energy, var- 
ious aspects have been investigated via the correlation 
functions, such as the dependences of the isospin of the 
emitting source 0, the impact parameter Q, the nu- 
clear symmetrical energy [l5l| and the total momentum 
of nucleon pairs and so on. The details about the 
EOS and the collision process can be revealed from the 
correlation functions well. 

Some groups have applied HBT technique to study 
the exotic nuclear reaction recently 0, 0, 0| . Stud- 
ies have been performed for many years for exotic nu- 
clei with the increasing availability of radioactive nuclear 
beams around the world. Among the various techniques 
to investigate on the exotic nuclei, the measurements 



of the total interaction cross section and the fragment 
momentum distribution of the projectile are the main 
methods to explore the exotic structure in the past years 
m El IBypjal mM il IliJglJsS. some nuclei 
like ^He^^/^^Li U, ^^^T^llllllJli are con- 
sidered as two-neutron halo ones and ^^Be H^jUll, ^^C 
[^Is^, ^^O [s^l are as one neutron halo, and the pro- 
ton halo structure has been also proved to exist in the 
structure of the nuclei ^B^,^'^N^, "^^Al [H etc. In 
terms of the structure of the halo, integral measurements, 
such as total reaction cross sections, are only sensitive to 
the overall size. Dissociation reactions, in which the core 
and/or nucleons are detected in the final state, can pro- 
vide some structure information (17ij-the major difficulty 
being the relationship between the initial and final states 
as dictated by the distorting effects of the reaction [l8l| . 
So it is very interesting to investigate the exotic nuclei 
via HBT technique further. As we know, the binding 
energy and the nucleon(s)-separation energy are impor- 
tant for the structure of the nuclei. The former indicates 
the stable level of the nucleus and the nucleon-nucleon 
relationship among the nucleus, and the later is a good 
criterion to verify the possibility of the exotic nucleus. 
These two have been studied through the density calcula- 
tion by RMF theory in the past years. Researches about 
these two factors via the nuclear collisions systematics are 
needed. In this Letter, we expect to explore the relation- 
ship between these two factors and the nucleon-nucleon 
correlation function value at very small relative momen- 
tum with help of Isospin-Dependent Quantum Molecular 
Dynamics (IDQMD) model which can describe the reac- 
tion dynamics on event by event basis. 



II. HBT TECHNIQUE 
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Firstly, we would like to recall the HBT technique. 
As we know the wave function of relative motion of 
light identical particles, when emitted in close proxim- 
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ity in space-time, is modified by the final-state interac- 
tion (FSI) and quantum statistical symmetries (QSS), 
and this is the principle of the intensity interferometry, 
i.e. HBT. The correlation function is defined as the ratio 
between the measured two-particle distribution and the 
product of the independent single-particle distributions: 



C{pi,p2) = 



dn^ / dpidp2 
dn / dpidn / dp2 



(1) 



where d^n/ dpidp2 represents the correlated two-particle 
distribution and dn/dpi and dn/dp2 is the independent 
single-particle distributions of particle 1 and 2, respec- 
tively. Usually, the projection C{q) = cqj^ onto the 

relative three-momentum [q = ^\pi — p2\) is used as the 
momentum correlation function, where the measured dis- 
tribution of pairs ( N(q) ) is divided by a reconstructed 
distribution of uncorrelated pairs ( D(q) ). cq is a normal- 
ized coefficient so that C{q) tends to 1 at high relative 
momentum, where the effects of FSI and QSS vanish. 
The deviation of C(g) from 1 thus reflects the informa- 
tion of the emission source. Other effects, arising from 
the form of the single-particle distributions or the exper- 
imental acceptances, are eliminated by the denominator 
ofEq. (1). 

Interpretation of correlation functions measured in 
heavy-ion collisions requires understanding the relation- 
ship of the parameters extracted from fitting the data and 
the true single-particle distributions at freeze-out. This 
relationship can be established by using an event gen- 
erator that models the collision dynamics, particle pro- 
duction, and then constructing a two-particle correlation 
function. The event-generator correlation functions are 
constructed from the positions and momenta represent- 
ing the single-particle emission distribution at the time 
of the last strong interaction, i.e. at freeze-out. The 
event-generator used in our task is the Isospin-dependent 
Quantum Molecular Dynamics transport model, which 
has been applied successfully to the studies of the heavy- 
ion collisions which includes the simulation of the reac- 
tion process 'SmiHlllElillil. For completeness, 
we would like to make a brief introduction to this model 
in following section. 



three components, the mean field, two-body collisions, 
and pauli blocking. Therefore, for an isospin-dependent 
reaction dynamics model it is important for these three 
components to include isospin degrees of freedom. What 
is more essential, in initialization of projectile and tar- 
get nuclei, the samples of neutrons and protons in phase 
space should be treated separately because there exists a 
large difference between neutron and proton density dis- 
tributions for nuclei far from the /3-stability line. Particu- 
larly, for neutron-rich nucleus one should sample a stable 
initialized nucleus with neutron-skin structure and there- 
fore one can directly explore the nuclear structure effects 
through a microscopic transport model. The IDQMD 
model has been improved based on the above ideas. 

In the IDQMD model, the nuclear mean field can be 
parametrized by 



C/(p,r,)=a(^) + /3(^r 



jPn " Pp) 
PO 



T, + U 



Yuk 



(2) 



with pq the normal nuclear matter density (here, 0.16 
a = -124 MeV, /? = 70.5 MeV and 7 = 2.0. p, 
Pn and Pp are the total, neutron, and proton densities, 
respectively; The coefficients a and /3 are parameters for 
nuclear equation of states (EOS), is the zth component 
of the isospin degree of freedom, which equals 1 or -1 
for neutrons or protons, respectively; Vc is the coulomb 
potential; and [7^"'= is Yukawa (surface) potential which 
has the following form: 



^Yuk ^ i:i —exp(Lm'^)\exp(-mn,)erf(VLr. 



)erf{^m + nj/V^)] (3) 

with Vy = 0.0074 GeV, m = 1.25 fm,-^ and L = 2.0 
frri^. The relative distance Vij — {fl — fj\. In our task, 
the so-called soft EOS with an incompressibility of K = 
200 MeV is used and the symmetry strength C = 32 MeV 



III. MODEL DESCRIPTION 

The Quantum Molecular Dynamics (QMD) approach 
is an n-body theory to describe heavy ion reactions from 
intermediate energy to 2 GeV/n. It includes five im- 
portant parts: the initialization of the target and the 
projectile; the propagation in the effective potential; the 
collisions between the nuclcons; the Pauli blocking ef- 
fect and the numerical tests. A g eneral review about the 
QMD model can be found in The IDQMD model is 
based on the QMD model affiliating the isospin factors. 

As we know, the dynamics in heavy-ion collisions 
(HIC) at intermediate energies is mainly governed by 



The NN cross section is the experimental parametriza- 
tion which is isospin dependent. Recently, studies of col- 
lective fiow in HIC at intermediate energies have revealed 
the reduction of the in- medium NN cross sections 
An empirical expression of the in-medium NN cross sec- 
tion [4JI is used: 



'NN 



^ — U + / )^NN 

Pa 



(4) 



0.2 which has been found to 

free 



with the factor / ; 

better reproduce the fiow data |43|. Here aj. 
the experimental NN cross section |43- The neutron- 
proton cross section is about three times larger than the 
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neutron-neutron or proton-proton cross section below 300 
MeV/nucleon. 

The Pauli blocking effect taken in IDQMD model is 
treated separately according to the neutron and the pro- 
ton: Whenever a collision has occurred, in the phase 
space we assume that each nucleon occupies a six- 
dimensional sphere with a volume of ?i'^/2 (considering 
the spin degree of freedom), and then calculate the phase 
volume, V, of the scattered nucleons being occupied by 
the rest nucleons with the same isospin as that of the 
scattered ones. We then compare 2Y /Tc' with a random 
number and decide whether the collision is blocked or 
not. 

For the initialization of the nucleons of the target and 
projectile, the IDQMD model distinguishes the proton 
and neutron from each other. The neutron and the pro- 
ton density distribution are determined from the Skyrme- 
Hartree-Fock (SHF) method with parameter set SKM* 
which can give reasonable density distribution for stable 
and neutron-rich nuclei |46| . One can obtain the radial 
positions of nucleons in the initial nuclei in terms of the 
Monte-Carlo method. In the model, the radial density 
can be written as: 



p{r) 



L 



2L ' 2rr. 
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FIG. 1: The solid circles with error bars are the two halo- neutron 
correlation functions (hh) in the collision of ^"^Be fragmented into 
12 Be -I- 2n at 35 MeV/n and the target is ^^cjig. The solid Une 
is the calculated two-halo neutrons correlation function. 



x[exp{^) - exp{ — j^)]. (5) 

And the momentum distribution of nucleons is gener- 
ated by means of the local Fermi gas approximation. The 
local Fermi momentum is given by: 

P^r) = h{37T^p,{r))^,{t = n,p). (6) 

The stability of the propagation of the initialized nu- 
clei has been checked in details and can last at least 200 
fm/c according to the evolutions of the average binding 
energies and the root mean square radii of the initialized 
nuclei. 

IV. RESULTS AND DISCUSSIONS 

Among the exotic nuclei, ^^Be nucleus, with four pro- 
tons and ten neutrons has received the most attentions, 
both theoretically and experimentally, due to its rather 
unique structure 0, 0; HE IHl • To test our approach, 
firstly we have analyzed the reaction of -'^^Be fragments 
into i^Be -t- 2n at 35 MeV/nucleon with the target ^^C. 
Fig. 1 shows the experimental and the calculated corre- 
lation functions for neutron and neutron. The two halo 
neutrons in calculation are defined as the emitted neu- 
trons in coincidence with ^^Be core. The solid line is 
the calculated two halo neutron correlation function. It 
shows clearly that the correlation function between the 
two halo neutrons reproduces the experimental data ex- 
cellently, which is consistent with the small two neutron 



separation energy. This directly leads to the neutron den- 
sity distribution of ^"^Be extending very large which is a 
major reason for the abnormal larger total reaction cross 
section and the ver y n arrow momentum distribution of 
the fragment ^^Be ^5^. In our work, the nucleons are 
defined as emitted if they do not belong to any clusters 
{A > 2) which are recognized by a simple coalescence 
model: ie. nucleons are considered to be part of a clus- 
ter if in the end at least one other nucleon is closer than 
rmin < 3.5 fm in coordinate space and Pmip. < 300 MeV/c 
in momentum space in the final state |^ |^ . This def- 
inition of nucleon emission is crucial for the good repro- 
duction of the experimental data in Fig. 1. The ^^_Be is 
therefore a good test case for the HBT method. 

Based upon the achievement of halo neutron - halo 
neutron correlation function of HBT results in IDQMD 
with the data, we further explore the proton-neutron cor- 
relation function as a function of binding energy or sepa- 
ration energy. The target is ^^C and the projectile is Be 
isotope. Only those events in which the neutron and the 
proton are emitted in the same event are accepted. The 
calculated results are shown in Fig. 2. The figure shows 
the proton-neutron correlation function for different Be 
isotopes and the insert of Fig. 2 shows the relationship 
between the strength of proton-neutron correlation func- 
tion CpN at 5 MeV/c and the binding energy per nucleon 
of the projectile Eunding- The solid line of the insert is 
just a linear fit to guide the eyes. 

One can notice that the behavior of the correlation 
functions between the proton and the neutron as a func- 
tion of neutron number and binding energy in Fig. 2. 
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FIG. 2: The proton-neutron correlation function for different Be 
isotopes. The meaning of symbols are illustrated in left bottom 
corner of Figure. The insert shows the relationship between the 
strength of proton- neutron correlation function Cpj^ at 5 MeV/c 
and the binding energy per nucleon of the projectile Ei^ij^^i^g. The 
collisions were simulated at 800 MeV/n of the incident energy and 
head-on collisions. The target is ^■^C. 



Generally the strength of HBT at very small relative mo- 
mentum shows a clear dependence on the neutron num- 
ber. The tendency of the CpM rises with the increasing 
Ebinding- As we know, among the projectiles we studied, 
the number of the protons is 4 and that of the neutron 
are gradually increasing, and this will reflect the stability 
of the nuclei. With the increasing of the mass number, 
the mean relationship between the nucleons will change 
weaker. Since the strength Cpjv of correlation function 
symbolizes the mean relationship between the emitted 
proton and neutron and the binding energy per nucleon 
associated with the tightness between the nucleon, the 
tendency shown in the insert of Fig. 2 reflects that the 
Cpjv can reveal the compactness of the nuclei. 

On the other hand, we also study the correlation func- 
tions between the proton and the most outside neutron. 
Here the most outside neutron is defined as the one which 
is the most far away from the spatial center at the FSI 
and it is also called as a halo neutron for simplification, 
even though it is not strict in physics sense. To obtain 
reasonable HBT results for halo neutron and proton, only 
those events in which the halo neutron and proton are 
emitted in the same event are accepted to investigate 
such a correlation function. The similar correlation func- 
tion was obtained as Cp^r and the relationship between 
the strength of proton-halo neutron correlation function 
CpH at 5 MeV/c and the single- neutron separation en- 
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FIG. 3: The relationship between the proton-halo neutron corre- 
lation function CpH at 5 MeV/c and the single-neutron separation 
energy of the nuclei. The solid line is just a linear fit. The collisions 
were simulated at 800 MeV/n of the incident energy and head-on 
collisions. The target is ^^C. 



ergy of the projectile Esep was extracted. The symbols 
of Fig. 3 show the calculated result and the solid line is 
just a linear fit to guide the eyes. It looks that, with the 
increasing of the Egep, CpH at 5 MeV/c rises gradually. 
Since the strength Cpu of correlation function between 
neutron and proton increases with the decreasing of the 
source size, the above Esep dependence of Cpu reflects 
that the emission source size of the neutron and proton 
shrinks with the separation energy of the single neutron, 
which is consistent with the extent of binding of single 
neutron via Eg^p- 

One should keep in mind that the Egep is thought as 
one of the important information for halo structure, for 
Be isotope, the nucleus ^^Be is considered to exist the 
halo configuration according to its extremely small one- 
neutron separation energy. Compared with other nuclei 
in Be isotopes, the Cpu at small relative momentum of 
it is smaller, i.e. there exists a larger spatial diffusion. 
That indicates that HBT of proton and halo neutron can 
show the different configuration of the nuclei Be isotopes 
well. These observations might illustrate that the CpH 
will symbolize the relationship between the most out- 
side neutron and the emitted proton, which is related 
to the extension level of the mass density distribution 
directly. Unfortunately, in the practical point of view, 
this correlation function Cpu is almost impossible to be 
measured. However, it still bear some information in the- 
oretical point of view. 



V. SUMMARY 

In summary, the intensity interferometry (HBT) tech- 
nique has been applied to investigate its sensitivity to 
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the binding energy and separation energy of neutron- 
rich nuclei from the break-up of nuclei by convoluting 
the phase-space distribution generated with the IDQMD 
model. Firstly we gave a well-fitted halo neutron - halo 
neutron correlation function from the break-up of ^'^Be 
on C target. Based upon this achievement of the good fit, 
we explore the dependence of the proton-neutron corre- 
lation function (Cpn) at small relative momentum with 
the binding energy {Ebmding) for Be isotopes. It was 
found that the correlation strength of Cpjv at small rel- 
ative momentum rises with the the binding energy. This 
changeable tendency of CpN with E^inding is here re- 
ported for the first time and it might be a potential good 
way to study the structure of the nuclei. Moreover, the 
proton- halo neutron correlation function {Cph) is also 
constructed from the break-up reactions. There exists 
the similar relationship between the Cpu at small rela- 
tive momentum and the separation energy (Esep) of Be 
isotopes as the relationship of Cpjv to E^inding- From 
theoretical point of view, CpH at small relative momen- 
tum is sensitive to Eg^p and this can be attributed to the 
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